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ABSTRACT 

After almost 2.5 years of actively accreting, the neutron star X-ray transient and eclipsing 
binary MXB 1659-29 returned to quiescence in 2001 September. We report on a Chandra obser- 
vation of this source taken a little over a month after this transition. The source was detected 
at an unabsorbed 0.5-10 keV flux of only (2.7 — 3.6) x 10 -13 erg cm -2 s _1 , which implies a 
0.5-10 keV X-ray luminosity of approximately (3.2 — 4.3) x 10 33 (d/10 kpc) 2 erg s _1 , with d 
the distance to the source in kpc. Its spectrum had a thermal shape and could be well fitted by 
either a blackbody with a temperature kT of ~ 0.3 keV or a neutron star atmosphere model with 
a kT of ~ 0.1 keV. The luminosity and spectral shape of MXB 1659-29 are very similar to those 
observed of the other neutron star X-ray transients when they are in their quiescent state. The 
source was variable during our observation, exhibiting a complete eclipse of the inner part of the 
system by the companion star. Dipping behavior was observed before the eclipse, likely due to 
obscuration by an extended feature in the outer part of a residual accretion disk. We discuss our 
observation in the context of the cooling neutron star model proposed to explain the quiescent 
properties of neutron star X-ray transients. 

Subject headings: accretion, accretion disks — stars: individual (MXB 1659-29) — X-rays: stars 



1. Introduction 

During outburst episodes, neutron star X-ray 
transients can be detected at luminosities of ~ 
10 36-38 erg s _1 (e.g., Chen, Shrader, & Livio 
1997). During those outbursts, the transients are 
very similar to the persistent sources with respect 
to their X-ray properties. The high X-ray luminos- 
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ity is very likely due to the accretion of matter onto 
the neutron star. These transients are character- 
ized by their bright outbursts, but most of the time 
they spend in a quiescent state in which they are 
orders of magnitude dimmer at all wavelengths. 
Fortunately, using sensitive imaging instruments, 
we are still able to detect them at X-ray luminosi- 
ties of ~ io 32 ~ 34 erg s _1 (e.g., van Paradijs et al. 
1987; Asai et al. 1996, 1998). The high sensitiv- 
ity camera's aboard Chandra and XMM-Newton 
are well suited to detect quiescent systems and 
obtain good X-ray spectra for the brightest sys- 
tems (see, e.g., Daigne et al. 2002; in 't Zand et 
al. 2001; Rutledge et al. 2001a, 2001b; Wijnands 
et al. 2001b, 2002b). To explain the low quies- 
cent X-ray properties, several models have been 
developed. For example, the X-rays could be due 
to the residual accretion of matter onto the neu- 
tron star or magnetospheric boundary, or the pul- 
sar emission mechanism might be active (see, e.g., 
Stella et al. 1994; Corbet 1996; Campana et al. 
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1998b; Menou et al. 1999; Campana & Stella 2000; 
Menou & McClintock 2001). Currently the most 
successful model is that in which the X-rays are 
due to the thermal emission from the neutron star 
surface, which will be referred to as 'the cooling 
neutron star model'. 

1.1. The cooling neutron star model 

In the cooling neutron star model (e.g., van 
Paradijs et al. 1987; Campana et al. 1998b; 
Brown, Bildstcn, & Rutledge 1998 and references 
therein) the emitted radiation below a few keV 
is thermal emission originating from the neutron 
star surface. Brown et al. (1998) argued that the 
neutron star core is heated by the nuclear reac- 
tions occurring deep in the crust when the star is 
accreting and this heat is released as thermal emis- 
sion during quiescence. If the quiescent emission is 
dominated by the thermal emission of the cooling 
neutron star, then the quiescent luminosity should 
depend on the time averaged (over 10 4-5 years) ac- 
cretion luminosity of the system (Campana et al. 
1998b; Brown et al. 1998). Thus, the quiescent lu- 
minosities of the detected systems can directly be 
compared with the predicted ones obtained from 
estimates of the long term accretion history of the 
sources. 

The neutron star cooling model also gives clear 
predictions for the spectral shape of the quiescent 
X-ray spectrum, which should be thermal. Al- 
though a simple blackbody model can be fitted to 
the data, the obtained radii of the emitting regions 
(of the order of only a few kilometers) are consid- 
erably lower than the predicted radii of neutron 
stars (Shapiro & Teukolsky 1983). To circumvent 
this discrepancy, it has been proposed that quies- 
cent neutron star systems do not emit a true black- 
body spectrum but a modified one (Brown et al. 
1998). When using blackbody models to fit such 
modified spectra, the effective temperatures will 
be overestimated and the emitting areas underes- 
timated. By fitting more realistic models, such 
as the so-called 'neutron star atmosphere models' 
(the non-magnetic models are appropriate for qui- 
escent neutron star systems; e.g., Zavlin, Pavlov, 
& Shibanov 1996 ), to the X-ray data, emitting 
radii were obtained which are consistent with the 
expected radii of neutron stars (Rutledge et al. 
1999, 2000). 

The cooling neutron star model cannot fully 



explain all characteristics of the quiescent emis- 
sion. For example, the power-law shaped spectral 
component which dominates the quiescent spectra 
above a few keV in several systems 7 (e.g., Asai et 
al. 1996, 1998; Campana et al. 1998a) cannot be 
explained by the cooling models. It is conceivable 
that this component might be described by one 
or more of the alternative models discussed above 
(in particular the residual accretion model) . How- 
ever, the observational results on this component 
and our understanding of its nature are very lim- 
ited. 

1.2. The quasi-persistent X-ray transients 

Recently, a sub-group of neutron star X-ray 
transients has received extra attention because 
of their potential to test the cooling neutron star 
model and to determine some of the physical prop- 
erties of the neutron star crust and core. These 
particular transients do not have traditional out- 
bursts which only last weeks to at most a few 
months, but instead they stay active for several 
years to over a decade (and maybe even longer). 
These systems have been called long-duration 
transients or quasi-persistent sources (e.g., Wij- 
nands et al. 2001b, Wijnands 2002). The long 
outburst behavior of those sources might be re- 
lated to the extended episodes (several months to 
several years) of low-level activity seen in other 
transients usually after they have exhibited bright 
outbursts (e.g., in 4U 1630-47, Aql X-l, 4U 1608- 
52, or SAX J1808.4-3658; Kuulkers et al. 1997; 
Bradt et al. 2000; Wijnands et al. 2001c; Wachter 
et al. 2002), although those episodes are generally 
less luminous (< 10 36 erg s _1 ) than the outbursts 
of the quasi-persistent sources (10 36 — 10 37 erg 

B" 1 ). 

In 'ordinary' (i.e., short-duration) transients, 
the accretion of matter will have only a very mi- 
nor effect on the thermal state of the crust, but for 
these quasi-persistent sources the prolonged accre- 
tion episodes can heat the crust to high temper- 
atures, considerably higher than that of the neu- 
tron star core (see Rutledge et al. 2002). When 
those systems become quiescent again, it might 



Note that not in all detected quiescent systems this power- 
law component could be detected and that the flux ratio 
of the power-law component with the thermal component 
varies considerably between sources. 
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take years to decades for the crust to return to 
thermal equilibrium with the core and the ini- 
tial quiescent properties of those systems might be 
dominated by the crust emission and not by the 
state of the core (as is the case in ordinary tran- 
sients). Monitoring observations of those systems 
in quiescence might even allow one to follow the 
cooling of the crust from which the heat conduc- 
tivity of the crust can be determined (Rutledge et 
al. 2002). 

Recently, one of the quasi-persistent systems 
(KS 1731-260) suddenly turned off after having 
actively accreted for over 12.5 years. A Chandra 
observation taken a few months after this transi- 
tion showed the source at a 0.5-10 keV luminos- 
ity of ~ 10 33 erg s" 1 (Wijnands et al. 2001b). 
An XMM-Newton observation of this system per- 
formed about half a year after the Chandra ob- 
servation, showed that the system had declined 
by a factor of ~3 in luminosity (Wijnands et al. 
2002b). If the quiescent emission from this sys- 
tem was dominated by the state of the crust, the 
decrease in luminosity within half a year strongly 
indicates that the crust must have a high heat con- 
ductivity (Wijnands et al. 2002b; using the cool- 
ing curves calculated for this system by Rutledge 
et al. 2002). In this scenario, the core temperature 
is expected to be lower than the crust temperature 
and the luminosity measured with XMM-Newton 
can be used as an upper limit on the core lumi- 
nosity. The quiescent luminosity of KS 1731-260 
is much lower than expected from its long term 
accretion history, and can only be explained in 
terms of the standard cooling model if this sys- 
tem is dormant for at least several thousand of 
years between outbursts (assuming all outbursts 
of this system are very similar to the last one, 
which might not be a valid assumption; Wijnands 
et al. 2001b; Rutledge et al. 2002). Alterna- 
tively, enhanced cooling processes might be active 
in the core, rapidly cooling it (e.g., Wijnands et 
al. 2001b, 2002b).' 

1.3. MXB 1659-29 

In 2001 September, the opportunity arose to 
use another quasi-persistent system to test the 
cooling neutron star model: MXB 1659-29. This 
source is an X-ray transient and was discovered in 
1976 by Lewin, Hoffman, & Doty (1976) during 
type-I X-ray bursts, which clearly demonstrates 
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Fig. 1.— The RXTE/ ASM light curve of MXB 
1659-29 clearly showing the 1999-2001 outburst. 
The time of our Chandra observation is indicated 
by the solid line. Note that our Chandra observa- 
tion was performed at times that the RXTE/ ASM 
(and also the RXTE/PCA) could not detect the 
source anymore (see text). 



that the compact object in this system is a neu- 
tron star. The source was detected several times 
between 1976 October and 1978 September with 
SAS3 and HEAO (Lewin et al. 1978; Share et al. 
1978; Griffiths et al. 1978; Cominsky, Ossman, & 
Lewin 1983; Cominsky & Wood 1984, 1989) and 
irregular X-ray variability was found in this system 
(Lewin 1979; Cominsky et al. 1983). Cominsky 
& Wood (1984, 1989) reported on the discovery 
of eclipses every ~7.1 hours, which can be identi- 
fied with the orbital period of the system. MXB 
1659-29 is one of only several non-pulsating neu- 
tron star low-mass X-ray binaries for which total 
eclipses have been observed (The other confirmed 
eclipsing neutron star systems are EXO 0748-676, 
GRS 1747-312, and AX J 1745. 6-2901, although 
several other systems [e.g., 4U 2129+47] show par- 
tial eclipses; Parmar et al. 1986; In 't Zand et 
al. 2000; Maeda et al. 1996; McClintock et al. 
1982). During later observations, using a variety 
of satellites (e.g., Hakucho, EXOSAT, ROSAT), 
the source could not be detected anymore (Comin- 
sky et al. 1983; Vcrbunt 2001). The pointed 
ROSAT observations in the early 1990's failed to 
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detected the source with a 0.5-10 keV upper limit 
on the unabsorbed flux of (1 — 2) x 10~ 14 erg cm~ 2 
s- 1 (Vcrbunt 2001; Wijnands 2002; Oosterbroek 
ct al. 2001). 

The source remained dormant until 1999 April, 
when in 't Zand et al. (1999) reported it to be ac- 
tive again in observations obtained with the Bep- 
poSAX Wide Field Camera. Wachter, Smale, & 
Bailyn (2000) and Oosterbroek ct al. (2001) ob- 
tained an updated ephemeris for the orbital pe- 
riod (using data obtained with the Rossi X-ray 
Timing Explorer [RXTE] and BeppoSAX). Stud- 
ies of its X-ray spectrum were performed using 
BeppoSAX (Oosterbroek et al. 2001) and XMM- 
Newton (Sidoli et al. 2001), and, using RXTE 
data, Wijnands, Strohmayer, & Franco (2001a) 
found ~567 Hz oscillations during X-ray bursts. 
Those oscillations are likely related to the neutron 
star spin frequency. The source remained bright 
for almost 2.5 years before it became dormant 
again in 2001 September (Wijnands et al. 2002a). 
Because of its long outburst duration, MXB 1659- 
29 may be classified as a quasi-persistent source. 
We had a Cycle 3 Chandra TOO proposal ap- 
proved to obtain a quiescent observation of the 
next quasi-persistent source that could turn off, 
within a month after the transition. As part of 
this proposal, MXB 1659-29 was observed on 2001 
October 15 using Chandra for ~ 19 ksec. Here we 
report on this observation. 

2. Observation, analysis, and results 

The RXTE All Sky Monitor (ASM) light curve 
of MXB 1659-29 is shown in Figure 1 (see also 
Wijnands et al. 2002a). As can clearly be seen, 
MXB 1659-29 stayed active during its last out- 
burst for almost 2.5 years before it could not be de- 
tected anymore with the RXTE/ ASM at the end 
of 2001 August. During more sensitive observa- 
tions using the proportional counter array (PCA) 
aboard RXTE, the source could still be detected 
until 2001 September 7, but it was undetectable in 
observations performed on 2001 September 14, 24, 
and 30 (Wijnands et al. 2002a; with upper limits 
on the flux of 0.5-1 mCrab; 2-60 keV). 

After the non-dctcction using the RXTE/PCA 
on 2001 September 14, we triggered our Cycle 3 
Chandra proposal and a Chandra/ ACISS observa- 
tion on MXB 1659-29 was performed on 15 Octo- 
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Fig. 2.— The CTwmdra/ACIS-S 0.5-8 keV im- 
age of the field of MXB 1659-29. The source is 
clearly detected together with two other weak X- 
ray sources. 



bcr 2001 between 16:31 and 22:24 UTC for a total 
of ~18.8 ksec of on-source time. No background 
flares occurred during the observation so all data 
could be used. We used the CIAO tools (version 
2.2.1) and the standard Chandra analysis threads 8 
to analyze the data. During our observation the 
ACIS-S3 CCD was used with a 1/4 sub-array (re- 
sulting in a frame time of 0.8 seconds). This con- 
figuration was used to reduce the possible pile-up 
problems in case the source had exceeded a flux 
level of ~ 10 -12 erg cm~ 2 s _1 . As demonstrated 
below, only about 2% pile-up occurred during our 
observation because of the relatively low flux of 
the source. 

2.1. X-ray image and light curve of MXB 
1659 29 

The obtained 0.5-8 keV image of the region 
around MXB 1659-29 is shown in Figure 2. We 
used the tool wavdetect to search the complete chip 
for point sources. The three sources visible in Fig- 
ure 2 are detected (their coordinates are listed in 
Tab. 1) together with two additional, rather dif- 
fuse sources, which could be truly diffuse sources 
or due to an elevated background count rate at 

8 Listcd at http://asc.harvard.edu 
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Fig. 3. — The optical I band image from Wachtcr 
ct al. (2000) during outburst together with the 
Chandra positions of MXB 1659-29 and CXOU 
J170205.9-295619, using circles with radii of 0.5". 
Clearly, the Chandra position of MXB 1659-29 is 
fully consistent with the optical one. 



those positions. In this paper, we are only inter- 
ested in the properties of MXB 1659-29 so we do 
not further investigate the significance of those po- 
tential sources. We used the optical I-band images 
of Wachter et al. (2000) taken during outburst 
to determine possible optical counterparts for the 
detected Chandra sources. We tied these optical 
images to the 2MASS J image available for this re- 
gion; the resulting image, including the Chandra 
positional circles, is shown in Figure 3. Clearly, 
the Chandra position of MXB 1659-29 is consis- 
tent with its optical position during outburst, with 
an offset of only ~ 0.15", well within the bore sight 
error (Aldcroft et al. 2000) and the errors on the 
optical position. The positions listed in Table 1 
are adjusted for this small off-set. The two extra 
Chandra sources do not have counterparts in the 
optical image. 

Figure 4 displays the Chandra/ ACIS-S light 
curve. It clearly shows an eclipse, which is very 
likely due to the obscuration of the inner X- 
ray emitting region by the companion star (cycle 
30823 using the linear ephemeris of Oosterbroek 
ct al. 2001). The eclipse duration (842±90 sec- 



onds) and the egress and ingress times (< 247 
and < 400 seconds, respectively) have been calcu- 
lated using the method outlined by Nowak et al. 
(2002). The obtained eclipse duration is consistent 
with that observed during outburst (<~ 900 sec- 
onds; Wachter et al. 2000). Eclipses in quiescence 
have also been observed for the neutron star tran- 
sient 4U 2129+47 (Nowak et al. 2002) and from 
the likely neutron star transient X5 in the globular 
cluster 47 Tuc (Hcinke ct al. 2001, 2002). Besides 
the eclipse, the Chandra/ ACIS-S light curve also 
shows clear dipping behavior a few kiloseconds be- 
fore the eclipse, which again is very similar to 4U 
2129+47 and X5 in 47 Tuc. Those dips might 
be due to the same process causing the dipping 
behavior of the source when it is still actively ac- 
creting: the dips might be due to obscuration of 
the central X-ray source by a large structure in the 
outer accretion disk, possible the impact point of 
the accretion stream from the companion star and 
the accretion disk (see also Nowak et al. 2002 and 
Heinke et al. 2001, 2002 for a discussion about the 
dips for 4U 2129+47 and X5, respectively). The 
decrease in count rate would then be due to an in- 
crease in the absorption column in front of the cen- 
tral X-ray source. This hypothesis can be tested 
by examining the spectrum of the source in and 
outside the dip (see § 2.3). Note that this hypoth- 
esis requires that residual disks are still present in 
quiescence without accretion occurring onto the 
neutron stars. Indeed such residual disks are ex- 
pected from the disk instability model to explain 
X-ray outbursts (see Lasota 2001 for a review) and 
independent proof for such disks came from opti- 
cal observations from quiescent X-ray transients 
(e.g., McClintock & Remillard 2000). 

2.2. The X-ray spectrum of the persistent 
emission 

All X-ray spectra were extracted using a circle 
with 5" in radius on the position of the source as 
determined using wavdetect and the background 
data were obtained by using an annulus on the 
source position with an inner radius of 7" and an 
outer one of 20". All obtained spectra were re- 
binned using the FTOOLS tool grppha into bins 
with a minimum of 15 counts per bin and we fitted 
the spectrum using XSPEC (version 11.1; Arnaud 
1996). We extracted spectra for the part of the 
observation during which the source was relatively 
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Fig. 4, - The Chandra/ ACIS-S 0.5-8 keV light 
curve of MXB 1659-29. Clearly visible are the 
eclipse and dipping behavior (about 6 ksec before 
the eclipse) . It is indicated in the figure which part 
of the data has been used to extract the persistent 
spectrum of MXB 1659-29 (see Fig. 5) and the 
dip spectrum. Each point represents 400 seconds 
of data. 



stable (called the 'persistent spectrum'; see Fig. 4; 
<~ 13.2 ksec of data) and for the part during the 
dips (called the 'dip spectrum'; Fig. 4; a total of 
~3.8 ksec of data). 

Many single component models fit the persis- 
tent spectrum well; however, the quiescent spec- 
tra of neutron star X-ray transients are most of- 
ten fit with a blackbody model or a neutron star 
atmosphere model. Therefore, we concentrate on 
those two models, using the model from Zavlin 
et al. 1996 (the non-magnetic case) as the atmo- 
sphere model. Theoretically it is expected that 
the emerging spectrum should resemble that as- 
sumed by the atmosphere model (e.g., Brown et 
al. 1998), however, observationally both the atmo- 
sphere model and the black body model produce 
equally satisfactory fits to the data. Because also 
our data do not allow for a rejection of the black- 
body model, and to allow comparison with previ- 
ous results on the quiescent spectra of neutron star 
X-ray transients, we also report the blackbody re- 
sults. In certain systems, a power-law tail above a 
few keV was found, and although such a power-law 




Black-body model 




Neutron star atmosphere model 



Fig. 5.— The Chandra/ ACIS-S spectrum of MXB 
1659-29 of the persistent part of the data. In the 
top panel, the solid line represents the best black- 
body fit to the data, in the bottom panel the best 
neutron star atmosphere model fit. 



component was not required by the data, we fitted 
the spectra with the above two models including 
a power-law component with a photon index of 1 
or 2 to obtain an upper limit on this component. 

The spectral results are listed in Table 2 and 
the persistent spectrum is shown in Figure 5. 
The column density Ah was left free in the fits 
and the obtained value is in the range of what 
has been observed previously for this source dur- 
ing outburst using BeppoSAX and XMM-Newton 
(0.13 - 0.35 x 10 22 cm- 2 ; Oostcrbrock et al. 2001; 
Sidoli et al. 2001). Both the blackbody and the 
atmosphere model fit the data well and the best-fit 
temperature was ~0.3 keV for the blackbody fits 
and ~0.1 keV for the atmosphere model 9 . When 



9 Using the obtained spectral parameters in this paragraph 
and the 0.8 frame time of our observation, we have used 
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assuming a distance of 10 kpc 10 , the radius for 
the emitting area was only a few kilometers using 
the blackbody model which is lower than the the- 
oretical expected radius for a neutron star. The 
atmosphere model gave radii that were consistent 
with those expected if the emission arose from the 
neutron star and the complete surface was radiat- 
ing, favoring the atmosphere model. For larger as- 
sumed distances the radii would increase, although 
the distance has to be unrealistically large (> 20 
kpc) for the radius obtained via the black-body 
model to become consistent with theoretical ex- 
pectations. The atmosphere model does not allow 
for the distance to become much larger than as- 
sumed because the obtained radius will quickly be 
inconsistent with theoretical expectations. 

The obtained unabsorbed 0.5-10 keV fluxes 
were between 3.0 and 3.6 xl0~ 13 erg cm~ 2 
s _1 , resulting in a luminosity of (3.6 — 4.3) x 
10 33 (d/10 kpc) 2 erg s _1 . When including a 
power-law component in the fit, it could not be de- 
tected significantly and its 0.5-10 keV flux could 
be constrained to be less than about 25%-35% 
(depending on photon index) of the 0.5-10 keV 
blackbody flux or that of the atmosphere compo- 
nent. 

2.3. The spectrum during the dip 

Due to the low statistics of the dip data, the 
source spectrum during the dip could not be con- 
strained by using only the dip data. However, by 
assuming that the decrease in count rate is only 
due to an increase in absorption by obscuring ma- 
terial (as has been found during outburst; e.g., 
Oosterbroek et al. 2001; Sidoli et al. 2001), we 
fitted the persistent spectrum and the dip spec- 
trum simultaneously using the same blackbody or 
atmosphere model but with different column den- 



PIMMS to estimate the degree of pile-up in our spectrum. 
We found that about 2% of the photons should be piled- 
up and the effects on the X-ray spectral parameters were 
only marginal. Therefore, we have not corrected the X-ray 
spectrum for this effect. 

The distance toward the source is not well known. Us- 
ing the luminosity of the X-ray bursts detected during the 
1999-2001 outburst, Oosterbroek ct al. (2001) obtained a 
range of 11-13 kpc but Muno et al. (2001) reported 10 kpc. 
In the remaining we will assume 10 kpc, but we note that 
it is rather uncertain. When appropriate we will discuss 
the effect of this uncertainty on the interpretation of our 
results. 
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Fig. 6. — The error ellipse of the fitted column 
density of the persistent emission part of the data 
(iVy 0rslstcnt ) versus that of the dip part (A^| lp ). 
The model used in XSPEC was the bbodyrad 
model. The plus represents the best fit value and 
the ellipses are for 1, 2, and 3a confidence levels 
(for two parameters). The straight line represents 
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sities. The spectral results are also listed in Ta- 
ble 2. The spectral parameters of the persistent 
emission were consistent with those obtained when 
only fitting the persistent emission spectrum, al- 
though the measured fluxes tend to be slightly 
lower. The column density of the dip spectrum 
is higher (by a factor of two) than that obtained 
for the persistent emission. Figure 6 shows the er- 
ror ellipse of the column density of the persistent 
emission spectrum versus that of the dip spectrum. 
Clearly, at high significance it is rejected that the 
column density of the dip spectrum is identical to 
that of the persistent emission. A systematic trend 
is clearly seen that the dip column density is al- 
ways larger than that of the persistent emission 
(a similar result is obtained when the atmosphere 
model is used instead of the blackbody model). 
Our results show that our data are consistent with 
a scenario in which the intrinsic source spectrum 
remains identical throughout the Chandra obser- 
vation, and that the decrease in count rate during 
the dip is due to an increase in absorption, but we 
cannot exclude other mechanisms responsible for 
the count rate decrease. For example, the data are 
also consistent with a decrease in luminosity with- 
out a change in spectral shape or with a decrease 
of the temperature. 
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The observation of a dip in the quiescent X-ray 
light curve of MXB 1659-29 might imply that ab- 
sorption by the outer disk also could occur, albeit 
less strong, during the orbital phases where no dips 
are observe, possibly effecting the spectral fits of 
the persistent emission. However, during the out- 
burst phase of the source type-I X-ray bursts were 
observed and nearly-coherent oscillations during 
those bursts were seen (Wijnands et al. 2001a), 
implying that the neutron star is observed directly. 
Very likely, this will also be true during quies- 
cence. Furthermore, the measured column density 
toward the source from the persistent quiescent 
emission is rather low and consistent with other 
estimates and measurements (Dickey & Lockman 
1990; Oosterbroek et al. 2001; Sidoli ct al. 2001) 
making additional absorption by the outer disk 
also unlikely. Therefore, we conclude that very 
likely we see the neutron star directly at orbital 
phases outside the dip and the eclipse, and our 
spectral results can directly be compared with 
the results obtained from other quiescent neutron 
star transients which have lower binary inclination 
than MXB 1659-29. 

3. Discussion 

We have presented a Chandra observation of 
MXB 1659-29 performed ~5 weeks after the last 
clear detection of the source with the RXTE/PCA 
(indicating that at that time the source was still 
actively accreting; Wijnands et al. 2002a). During 
our Chandra observation, we detected the source 
at a luminosity of (3.2 - 4.3) x 10 33 (d/10 kpc) 2 erg 
s _1 , and its spectrum could be well described by a 
thermal component (either a blackbody model or 
a neutron star atmosphere model). The obtained 
luminosity and the shape of the X-ray spectrum of 
MXB 1659-29 resemble those obtained for other 
quiescent neutron star systems, strongly suggest- 
ing that the source was quiescent during our Chan- 
dra observation. 

3.1. The cooling neutron star model 

As argued by Rutledge et al. (2002), for sys- 
tems which are actively accreting for long periods, 
the crust might have been heated to very high tem- 
peratures (probably considerably higher than that 
of the core) and the quiescent properties might be 
dominated by the thermal state of the crust rather 



than that of the core. Although the duration of the 
accretion episode of MXB 1659-29 is considerably 
shorter (factor of -5) than that of KS 1731-260 or 
X 1732-304 (both of which had accretion episodes 
of more than a decade; Wijnands et al. 2001b; 
Wijnands et al. 2002b; see also Guainazzi, Par- 
mar, & Oosterbroek 1999 for X 1732 304), here 
we assume that the 2.5 year accretion episode of 
MXB 1659-29 has had a considerable effect on the 
state of the crust, similar to what has been argued 
for KS 1731-260 (Rutledge et al. 2002; although 
likely less extreme), . 

This assumption is also supported by the fact 
that in the early 1990's, the source could not 
be detected in quiescence using a ROSAT ob- 
servation and only a 0.5-10 keV upper limit of 
(1 - 2) x 10 32 (d/10 kpc) 2 erg s" 1 could be ob- 
tained (Verbunt 2001; Wijnands 2002; Ooster- 
broek et al. 2001) which is about an order of 
magnitude lower than the luminosity we have de- 
tected during our Chandra observation. There- 
fore, we can conclude that the Chandra quiescent 
luminosity is not the rock bottom quiescent lumi- 
nosity of this system and the quiescent properties 
for MXB 1659-29 during the Chandra observation 
were dominated by that of the crust and not by 
the core. As an upper limit on the flux due to 
the cooling neutron star core, we will assume the 
upper limit provided by ROSAT and this upper 
limit will be used to test the cooling neutron star 
model. 

In order to test the model, the time averaged ac- 
cretion rate has to be estimated. The last outburst 
was fully covered with the RXTEj ASM instru- 
ment (Fig. 1) and lasted for ~ 2.5 years. The 2-10 
keV luminosity during this state as obtained with 
BeppoSAX and XMM-Newton was about 6 x 10~ 10 
erg cm~ 2 s _1 (Oosterbroek et al. 2001; Sidoli et al. 
2001). Using the spectral model and the spectral 
parameters given by Oosterbroek et al. (2001), we 
inferred (by simulating the spectrum in XSPEC) 
that the bolomctric luminosity can be at least a 
factor of two larger. However, the luminosity it- 
self was variable (by a factor of a few; Fig. 1) dur- 
ing the outburst, and large uncertainties might be 
present in the 'typical' outburst bolometric lumi- 
nosity. However, in the rest of the discussion we 
assume that bolometric outburst flux was typically 
(5 — 10) x 10~ 10 erg cm~ 2 s _1 during the last out- 
burst. 
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The 1999-2001 outburst could be an atypi- 
cal one for MXB 1659-29 and previous outbursts 
might have been less bright and/or less long (i.e., 
more like those of the ordinary transients). We 
have searched the literature for reports on detec- 
tions of this source in the past and we found that 
the source was conclusively detected in X-rays in 
1976 October, 1977 June, July 11 , and September, 
and 1978 March and September using SAS3 and 
EE AO (Lewin et al. 1976; Lewin et al. 1978; 
Share et al. 1978; Griffiths et al. 1978; Lewin 
1979; Cominsky et al. 1983; Cominsky & Wood 
1984, 1989) and in the optical on 1978 June 1 and 
1979 June 27 - July 2 (Doxsey et al. 1979; Comin- 
sky et al. 1983). No information is available dur- 
ing the periods in-between those observations. Al- 
though it cannot be excluded that at those occa- 
sions, the source was in quiescence, we consider it 
unlikely that the source would only be active dur- 
ing times when it was observed with an X-ray or 
optical instrument and dormant when no instru- 
ment looked at the source. Therefore, it is likely 
that during the complete period from 1976 Octo- 
ber until early July 1979 the source was actively 
accreting for over 2.5 years, especially because the 
recent outburst had a similar duration. If true, 
then this would constitute the first indication that 
different outbursts of quasi-persistent sources may 
have similar durations and that the long duration 
of those outbursts might be a common property of 
those sources. 

The first reported non-detection of the source 
was on 1979 July 17-25 (Cominsky et al. 1983) in 
optical (V>22-23) and with Eakucho (no X-ray 
upper limits were provided). Prior to 1976, the 
source might have also been detected during the 
period 1971 to 1973 using Uhuru (classified as 4U 
1704-30; Forman et al. 1978), although this iden- 
tification with MXB 1659-29 is not certain and 
we will assume that they are two different objects 
(if 4U 1704-30 can be identified with MXB 1659- 
29 then the source might have exhibited an ex- 
tra outburst during that period or it might have 
been active for a period of ~7 years). The exact 

During 1977 June and July, the observations were not very 
sensitive to persistent X-ray emission and none was de- 
tected (e.g., Cominsky & Wood 1984, 1989). But X-ray 
bursts were observed from the source, indicating that also 
during those observations, the source was accreting, albeit 
at a low level. 



fluxes during the observations in the period Oc- 
tober 1976 to early July 1979 are in the range 1 
to 6 xl0~ 10 erg cm~ 2 s _1 , but they have large 
uncertainties because the exact energy range was 
not always quoted (if quoted it was 1-10 keV or 
2-10 keV), it was unclear if the fluxes were ab- 
sorbed or unabsorbed, and the assumed spectral 
shape was not always similar (often assumed to be 
Crab like) and quite different than observed with 
BeppoSAX and XMM-Newton during the 1999- 
2001 outburst (Oosterbroek et al. 2001; Sidoli et 
al. 2001). However, we will assume that all fluxes 
are for the 2-10 keV range and unabsorbed, and 
a bolometric flux of about twice the quoted values 
(as inferred above for the BeppoSAX results on 
MXB 1659-29). Therefore, during the outburst in 
the late 1970's, the source was actively accreting 
for a period of at least ~2.5 years at a bolometric 
flux level of 2 to 12 xlO 10 erg cm~ 2 s _1 . Although 
quite uncertain, this is remarkably similar to the 
values of the 1999-2001 outburst and for simplic- 
ity we assume that the typical outburst duration 
is 2.5 years and that the bolometric fluxes during 
outburst is 5-10 xlO -10 erg cm~ 2 s _1 . 

Using the Brown et al. (1998) model (assuming 
standard cooling processes), the predicted quies- 
cent flux F q for this source would then be (Wij- 
nands et al. 2001b; see also Rutledgc et al. 2002) 
F q w t ^ x -y^, with (F ) the average flux dur- 
ing outburst (5-10 xl0~ 10 erg cm~ 2 s _1 ), t a the 
average time the source is in outburst (2.5 years), 
and t q the average time the source is in quiescence 
(~21 year). This results in a predicted quiescent 
flux of 4-8 x 10 -13 erg cm -2 s _1 . Remarkably, this 
value is very similar to the quiescent flux observed 
during our Chandra observations. However, as ex- 
plained above, based on the ROSAT non-detection 
of the source, the core flux is likely at least an 
order of magnitude lower than this, which would 
make the predicted core flux considerably higher 
than that truly originating from the core. Similar 
to KS 1731-260 (Wijnands et al. 2001b, 2002b) 
this low core flux (and thus temperature) might 
be due to enhanced core cooling instead of the as- 
sumed standard core cooling in the Brown et al. 
(1998) model. 

Despite the fact that the last two outburst 
episodes are likely of similar duration, it cannot be 
excluded that they are not typical for the source 
and that most of the time MXB 1659-29 exhibits 
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short duration outbursts. If the typical outburst 
duration of MXB 1659-29 is not 2.5 years but in- 
stead 0.25 years (3 months) or shorter, then the 
predicted quiescent flux will be consistent (within 
the uncertainties of the model and assumptions) 
with the ROSAT upper limit. The higher Chan- 
dra quiescent luminosity is again due to the state 
of the crust which should be considerably heated 
during the long accretion episode. 

3.2. Crust cooling 

If during the Chandra observation the X-ray 
emission was dominated by thermal emission from 
the crust, then further quiescent observations of 
MXB 1659-29 will enable studies of the cooling of 
the crust in this system. The ROSA T flux upper 
limit suggests that the crust flux should eventually 
decrease to at least this level. When the crust will 
be thermally relaxed with the core, no significant 
further decrease of the quiescent flux is expected 
and from this bottom flux level the state of the 
core can be inferred from which the cooling mod- 
els can be better constrained. For KS 1731-260 
it had already been found that its quiescent lu- 
minosity decreased by a factor of 3 within half a 
year time (likely due to a temperature decrease), 
indicating a highly conductive crust (Wijnands et 
al. 2002b) . It would be of interest to determine if 
such a rapid cooling will also be observed for MXB 
1659-29 or if the neutron star crust in this system 
has a significantly lower conductivity. 

In the latter case, the system should be at the 
Chandra quiescent luminosity for several years to 
decades. Using the fact that 13 years after the 
1976-1978 outburst ROSAT observed an order of 
magnitude lower quiescent luminosity we might 
set already an upper limit on the crust cooling 
time. When assuming that shortly after the end 
of the 1976-1978 outburst the quiescent luminos- 
ity was similar to our measured Chandra luminos- 
ity the cooling time of the crust is at least about 
a factor of 10 in luminosity per decade. Due to 
differences in quiescent times, outburst times, and 
the time averaged accretion rates between MXB 
1659-29 and KS 1731-260 the cooling curves cal- 
culated for KS 1731-260 by Rutledge et al. (2002) 
cannot be used for MXB 1659-29. However, if 
those MXB 1659-29 cooling curves resemble those 
of KS 1731-260, then tentatively it might be con- 
cluded that also the neutron star crust in MXB 



1659-29 has a high conductivity (and enhanced 
core cooling is suggested). We await specifically 
calculated cooling curves for MXB 1659-29 and 
further monitoring observations using Chandra or 
XMM-Newton in order to be conclusive about the 
properties of the neutron star in this system. 

3.3. Very low thermal emission in quies- 
cent neutron star systems? 

For most quiescent neutron star X-ray tran- 
sients it has been inferred that the thermal emis- 
sion is in the range of a few times 10 32-33 erg s _1 . 
However, recently, indications have been found 
(using XMM-Newton data) that the thermal emis- 
sion from the accretion driven millisecond X-ray 
pulsar SAX J1808. 4-3658 in its quiescent state 
might be as low as a few times 10 30 erg s _1 (Cam- 
pana et al. 2002; who suggested enhanced core 
cooling for this low thermal luminosity) . Although 
the statistics of those results were not overwhelm- 
ing and have to be confirmed with additional ob- 
servations, it is an interesting possibility and if 
the processes which produce the quiescent emis- 
sion (which has a power-law shape; Campana et 
al. 2002) for this system would become inactive, 
SAX J1808. 4-3658 would become rather dim in 
quiescence. Such weak quiescent neutron star sys- 
tems might also be suggested by the indications 
found for enhanced cooling processes in certain 
systems, such as MXB 1659-29. The possibility 
of dim quiescent neutron star systems raises the 
question of how dim certain neutron star systems 
can become? 

The answers to this question will have impli- 
cations for our understanding of quiescent X-ray 
binaries. It has been found that those X-ray bina- 
ries which harbor a black hole instead of a neutron 
star can be at least one to two orders of magnitude 
less luminous in quiescence than the average neu- 
tron star system. This difference has been used as 
evidence that the black holes have event horizons 
(e.g., Garcia et al. 2001 and references therein) 
in contrast to the surfaces of neutron stars. How- 
ever, if certain neutron star systems might become 
similarly dim (e.g., due to enhanced core cooling 
or very dim outbursts), than this luminosity dif- 
ference will disappear and with it the evidence for 
event horizons in black hole systems. The quies- 
cent spectrum of SAX J1808. 4-3658 (Campana et 
al. 2002) indicates that the spectrum of such sys- 



10 



terns might not be dominated by a thermal com- 
ponent but might have a power-law shape, similar 
to what has seen for the black hole systems (Kong 
et al. 2002) removing also the spectral differences. 

Another area in which it might be important to 
determine the full range of the luminosity distri- 
bution of quiescent neutron star systems, is that 
of the study of the low-luminosity X-ray sources 
in globular clusters. Based on their luminosities 
and their X-ray spectra (when enough statistics 
are available), those sources which have luminosi- 
ties above a few times 10 32 erg s -1 , have been clas- 
sified as quiescent neutron star systems and those 
which have luminosities below 10 32 erg s^ 1 , as 
another type of object (possible cataclysmic vari- 
ables or millisecond radio pulsars). However, the 
observed low luminosity of 5 x 10 31 erg s" 1 of SAX 
J1808.4-3658 (Dotani, Asai, & Wijnands 2000; 
Campana et al. 2002) already shows (irrespective 
of what exactly causes the X-rays in this source) 
that those dim sources might be quiescent neutron 
star transients. The fact that the spectrum of SAX 
J1808. 4-3658 appears to be considerably harder 
(Campana et al. 2002) than the average quies- 
cent neutron star spectrum indicates that classifi- 
cations on hardness ratio might lead to erroneous 
results. The possibility that the thermal emis- 
sion of SAX J1808. 4-3658 might be very low even 
suggests that those sources which have luminosi- 
ties down to only a few times 10 30 erg s _1 might 
also be quiescent neutron star systems. Therefore, 
we conclude that classifying low-luminosity glob- 
ular cluster sources based on their luminosity and 
broad band spectral shape (i.e., hardness ratio's) 
might possibly lead to misleading results. The full 
extent of those errors depends on the full lumi- 
nosity distribution of the quiescent neutron star 
transients and how the spectrum correlates with 
the luminosity. Those properties have to be deter- 
mined before a complete picture of the nature of 
the low-luminosity globular cluster sources can be 
understood 12 . 
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Although we have only discussed quiescent neutron star 
systems, similar uncertainties in the luminosity distribu- 
tion of CVs or millisecond radio pulsars exist. Those un- 
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resolved before a full understanding of the low-luminosity 
globular cluster source population can be obtained. 
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Table 1 

Measured Chandra positions of the detected sources 



Name 






Coordinates 

R.A. 


Dec. 




MXB 1659-29 
CXOU J170207.1 
CXOU J170205.9- 


-295535 
-295619 


yjh 02 m 

yjh Q 2 m 
yjh Q 2 m 


06.54 s ±0.02 s -29° 
07.06 s ±0.02 s -29° 
05.93 s ±0.02 s -29° 


56' 44.1" 
55' 34.7" 
56' 19.2" 


±0.3" 
±0.3" 
±0.4" 



Table 2 
Spectral results 



Parameter 




Persistent 


Persistent + dip 




Blackbody 


Hydrogen atmosphere 


Blackbody 


Hydrogen atmosphere 


^persistent ^ q22 cm _ 2) 

< ip (10 22 cm" 2 ) 
kT (keV) 

Radius ( 10 d kpc km) 

■^persistent 

X 2 /dof 


22+ - 07 

u - zz -0.06 

0.28±0.02 
2 4+ 1 ' 5 

3.0 
32.2/37 


0-28+^ 

o-io±8:S 
i7if 

3.6 
31.0/36 


U - iy -0.06 

S5+ 013 

u - oo -0.08 
n O o+0.02 

2 1+ 1 - 6 

2.7 
45.2/43 


0.27±0.07 

n 44+°- 13 

u - li -0.05 

14t 26 
3.4 
42.6/42 



Note: The error bars represent 90% confidence levels. The fluxes are unabsorbed, in the 0.5-10 
keV range, and in units of 10" 13 erg cm" 2 s _1 . As blackbody model, we used the bbodyrad model 
in XSPEC. For the hydrogen atmosphere model, we used the model of Zavlin et al. 1996 (the non- 
magnetic case) and for this model the neutron star mass was fixed to 1.4 M and the temperature and 
radius arc for an observer at infinity. For both models the distance d was assumed to be 10 kpc. 
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